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FOREWORD 


This  report  describes  an  in-house  effort  conducted  by  personnel  of  the 
Rockwell  International  Science  Center,  under  Project  7351  (DARPA  order  2828)  of 
General  Order  595,  "Interdisciplinary  Program  for  Quantitative  Flaw  Definition." 

1  he  work  reported  herein  was  performed  during  the  period  1  April  1978 
to  31  August  1978,  under  the  direction  of  the  author,  Dr.  C.  V.  Vasile.  The 
report  was  released  by  the  author  in  November  1978. 

This  report  is  Part  I  of  two  parts  concerned  with  detection  of  cracks 
under  fasteners  with  EMATS  (electromagnetic  acoustic  transducers) . 
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I.  INTRODUCTION  AND  SUMMARY 


This  report  describes  a  five-month  stud;/  having  as  an  objective  the 
demonstration  of  feasibility  and  generation  of  design  data  for  an  ultrasonic 
technique  for  the  nondestructive  evaluation  (NDE)  of  fastener  holes  In  the 
lower*  Inaccessible  section  of  an  aircraft  wing  lap  joint.  The  resulting 
measurement  technique  utilizes  the  excitation  of  horizontally  polarized  shear 
(SH)  waves  at  250  KHz  (0.5  In.  (1.27  cm)  wavelength  In  aluminum)  with  non- 
contact,  periodic  permanent  magnet  electromagnetic  transducers  (EMATs)  and  a 
dedicated  electronics  test  set.  The  transducers  are  placed  on  the  accessible 
portion  of  the  wing  adjacent  to  the  fastener  hole  in  question  and  the  flaw 
Information  Is  separated  from  signals  Introduced  by  geometrical  complexities 
of  the  structure  by  real  time,  Fourier  analysis  techniques. 

More  specifically  a  pair  of  EMATs  are  placed  as  close  as  possible  to 
a  single  fastener  hole,  thereby  forming  an  "ultrasonic  filter"  structure.  The 
presence  of  the  hole  perturbs  the  passband  of  this  filter,  while  the  presence 
of  a  crack  Introduces  a  further  significant  change  In  both  the  amplitude  and 
phase  response  of  the  "filter."  In  order  to  make  this  observation  It  Is 
necessary  to  perform  a  real-time  Fourier  analysis  of  the  ultrasonic  signals. 
Accomplishing  this  task  required  considerable  effort  In  modifying  and 
Improving  the  performance  of  an  existing  correlation  receiver  and  dedicated 
test  set.  While  such  effort  detracted  from  the  amount  of  test  data  we  would 
otherwise  have  collected,  development  of  this  improved  electronic  signal 
processing  was  essential  In  gathering  meaningful  and  reproducible  data. 

In  summary,  we  have:  (1)  developed  a  real-time  Fourier  transform 
signal  processor,  (2)  used  It  to  gather  experimental  data  demonstrating  the 
feasibility  of  detecting  cracks  growing  from  fasteners  in  lap  joints  of  fixed 
geometry,  and  (3)  deduced  a  simple  ultrasonics  model  capable  of  explaining 
many  of  the  experimental  data.  The  primary  problem  which  remains  Is  the  fact 
that  changes  In  the  geometry  of  the  wing  lap  joint  produce  interfering 
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effects,  and  further  work  must  be  done  to  learn  to  generate  reliable  crack 
information  in  regions  where  such  changes  are  occurring  or  the  geometry  is 
unknown.  Based  on  these  results,  we  recommend  a  second  phase,  the  goal  of 
which  would  be  to  design  and  construct  an  Improved  test  Inspection  system,  use 
It  to  gather  a  more  extensive  baseline  data  set  on  a  statistically  meaningful 
number  of  samples,  and  analyze  the  data  to  develop  specific  approaches  for 
obtaining  flaw  Information  In  regions  of  geometrical  changes  and  to  verify  the 
reliability  of  the  resulting  algorithms. 
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II.  BACKGROUND 

The  problem  is  summarized  in  Fig.  1.  The  upper  half  of  a  lap  joint 
can  be  inspected  by  available  techniques,  such  as  the  Rotoscanner.  However, 
the  lower  half  is  not  as  easily  inspected.  This  has  been  attributed  to 
variations  in  the  ultrasonic  transmission  properties  of  the  thin  sealant  layer 
(a  matter  that  is  discussed  further  in  the  Appendix).  Our  approach  combines 
four  innovative  Ideas.  Noncontact  EMATs  are  used  to  excite  horizontally 
polarized  shear  waves  on  member  B.  These  propagate  around  the  corner  and 
enter  the  lap  joint  region.  There  they  are  scattered  by  the  hole,  crack 
combination,  and  this  information  propagates  back  around  the  corner  to  a 
receiving  transducer,  also  on  member  B.  Finally,  the  flaw  Information  Is 
deduced  from  the  other  geometrical  Information  by  signal  processing,  in  this 
case  near  real  time  Fourier  analysis  using  a  correlation  receiver.  The 
remainder  of  this  section  provides  background  data  which  illustrates  these 
ideas,  as  they  existed  at  the  time  of  Initiation  of  the  program. 

Fig.  2  illustrates  the  basic  principles  of  an  FMAT  (electromagnetic- 
acoustic  transducer),  a  device  which  can  excite  ultrasonic  waves  with 
no  couplant,  through  a  loosely  bonded  paint  layer.  In  its  most  elemen¬ 
tary  form,  it  consists  of  a  single  wire,  next  to  a  metal  part,  which 
excites  eddy  currents  at  the  desired  ultr, sonic  frequency.  If  a  static 
magnetic  field  is  also  provided,  these  eddy  currents  will  experience  Lorentz 
forces  (the  same  forces  which  drive  an  electric  motor),  and  ultrasonic  waves 
are  radiated.  These  waves  can  be  detected  by  reciprocal  processes, 

A  real  transducer  Is  composed  of  an  array  of  wires  and  magnets  which 
produce  the  forces  required  to  form  the  desired  beam.  One  such  configuration, 
known  as  the  PPM  (periodic  permanent  magnet)  is  shown  in  the  right  half  of 
Fig.  2.  This  excites  horizontally  polarized  shear  (SH)  waves  which  are 
polarized  parallel  to  the  metal  surface.  Such  waves  are  particularly  useful 
for  the  present  problem  because  they  do  not  mode  convert  upon  reflection  from 
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Detection  of  Cracks  under  Fasteners 


Fig.  2  Coup i ant  Free  Inspection  with  EMATs. 


parallel  surfaces.  They  can  only  he  excited  hy  EMATs,  since  fluid  coupled 
angle  beam  systems  dn  not  support  this  component  of  shear.  In  this  program, 
the  transducer  structure  has  evolved  to  one  further  level  of  sophistication, 
as  shown  in  Fig.  3.  Here  the  magnets  have  been  separated  into  two  oppositely 
polarized  stacks,  and  the  coil  is  In  the  form  of  a  hairpin  loop.  Since  both 
the  magnetic  bias  and  current  change  direction,  the  force  is  the  same  as  in 
the  previous  configuration.  However,  the  electrical  properties  of  the 
structure  are  greatly  improved. 

Fig.  4  illustrates  the  principles  of  guided  wave  propagation  wh;ch 
are  being  used.  There  are  a  number  of  ultrasonic  modes  which  can  propagate 
a'ong  a  flat  plate.  These  are  characterized  by  a  particular  profile  of 
displacement  and  stress  (lower  right)  and  whose  frequency  (f  =  u/Zit)  and 
wavelength  ( \ )  are  related  by  a  set  of  dispersion  curves  (lower  left).  These 
inodes  can  propagjte  around  corners.with  some  reflection,  and  thereby  sense 
flaws  in  the  lap  joint  region. 

Fig.  5  presents  a  set  of  experiments  which  were  performed  to  demon¬ 
strate  the  ability  to  detect  cracks  growing  from  holes  by  measuring  the 
scattering  properties  of  these  waves.  For  these  measurements,  the 
noridi spersi ve,  n  -  0,  mode  was  selected.  The  upper  two  sets  of  data  show  the 
scattering  measured  for  a  1/4  in.  (0.635  cm)  hole  in  the  center  of  a  1/8  in. 
(0.318  cm)  plate,  as  measured  at  500  KHz.  To  the  left  we  see  the  scattering 
properties  of  the  hole.  To  the  right  we  see  the  changes  produced  when  an 
0.3  in.  (0.768  cm)  saw  slot  was  placed  in  the  hole.  The  bottom  two  sets  of 
data  show  the  results  that  were  obtained  when  the  hole  was  near  the  edge, 
rather  than  in  the  center,  of  the  plate.  The  extra  reflection  from  the  edge 
complicates  the  signals,  but  again  the  cracked  hole  is  clearly 
distinguishable. 

In  order  to  separate  the  flaw  Information  from  that  due  to  other 

geometrical  effects,  signal  processing  is  required.  Fig.  6  illustrates 
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Fig.  5  Scattering  from  Cracked  Hole. 
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several  candidate  procedures.  In  the  upper  left  Is  shown  a  plot  of  the 
nagnitude  of  the  backscattered  signal  from  a  slotted  hole.  The  phase 
information,  suppressed  in  that  presentation,  is  revealed  in  the  display  to 
the  right.  Here, both  the  In-phase  (I)  and  out-of-phase  (Q)  components  of  the 
signal  are  shown  as  a  function  of  angle.  These  may  be  thought  of  as  the  sine 
and  cosine  transforms  of  the  signal.  A  possibly  more  useful  representation  of 
this  information  Is  presented  at  the  bottom  left.  Here  the  I  component  Is 
plotted  versus  the  Q  component  with  angular  orientation  as  a  parameter.  The 
resulting  display  Is  quite  similar  to  that  which  has  proved  so  powerful  In  the 
presentation  and  Interpretation  of  eddy  current  data. 

The  remainder  of  the  report  describes  how  these  elements,  EMATs,  SH 
modes  propagating  around  corners,  scattering,  and  signal  processing,  have  been 
combined  to  develop  a  candidate  solution  to  the  problem. 
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III.  SYSTEM  CONCEPT 

Fig.  7  illustrates  a  representative  cross  section  of  the  structure 

to  be  inspected.  The  troublesome  area  Is  the  lower  lip  of  the  joint  as  It 
cannot  be  accessed  from  below  in  practice.  The  approach  we  selected  was  to 
locate  a  transmitting  EMAT  on  member  B  and  direct  ultrasonic  energy  In  the 
form  of  the  lowest  order,  horizontally  polarized  plate  mode  (e.g.,  the 
n  »  0  SH  mode)  past  the  fastener  and  measure  the  phase  and  amplitude 
characteristics  of  the  transmitted  signal  as  detected  by  a  similar  receiving 
EMAT. 


The  dimensions  illustrated  In  Fig.  7  are  merely  representative  of 
actual  structures  and  were  used  In  the  fabrication  of  a  mock-up  measuring 
24  x  24  in.  (61  x  61  cm)  containing  19  D  size  drill  holes  (.246  in.  (0.62  cm) 
dla)  spaced  1-1/4  In.  (3.18  cm)  on  centers.  Actual  wing  sections  have  tapered 
fastener  holes,  nonuni  form  hole  spacings  as  well  as  nonuniform  dimensions 
and  sealant  in  the  joint.  Therefore,  it  Is  Important  to  design  a  test  system 
that  can  minimize  and/or  account  for  these  variable  geometrical  parameters  In 
order  to  avoid  obscuring  the  presence  of  a  crack  which  is  the  geometrical 
discontinuity  which  we  are  seeking. 

Fig.  8  illustrates  the  location  of  a  pair  of  EMATs  and  some  impor¬ 
tant  path  lengths  In  the  configuration  that  was  ultimately  selected.  Note 
the  symmetrical  location  of  transmitter  and  receiver  and  the  fact  that  a 
narrow  wing  section  can  be  used  since  the  back  radiation  patterns  of  the  EMATs 
are  not  directed  towards  each  other,  even  In  the  presence  of  a  back  reflecting 
edge.  This  Is  Important  for  practical  reasons.  As  Illustrated,  Identical 
EMATs  (with  the  exception  of  different  matching  transformers)  are  used  for 
transmitter  and  receiver.  The  three  ultrasonic  wave  paths  Indicated  In  Fig.  2 
are  drawn  from  the  center  of  phase  of  each  four-period  long  transducer  and  are 
key  In  the  simple  model  to  be  described  below. 
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Fig.  7  Representative  Cross-Section  of  Wing  Lap  Joint 


Fig.  8  Location  of  SH  Wave  EMATs  with  Respect  to  Fastener  Hole. 


The  transducer  configuration  was  arrived  at  after  numerous 
theoretical  and  experimental  Investigations.  Early  In  the  program  we  used 
transducers  having  a  .25  In.  (0.64  cm)  period  which  therefore  operated  on  a 
9/32  In.  (0.71  cm)  thick  aluminum  plate  at. 488,  .54,  .645,  .78,  .92,  and 
1.1  MHz  for  the  n  ■  0,  1,  2,  3,  4,  5  SH  modes,  respectively.  Early  In  the 
program  a  large  wing  section  mock-up  was  not  available  and  a  broadband  system 
was  therefore  set  up  In  order  to  make  measurements  on  flat  plates,  right  angle 
bends,  and  jogs  with  all  of  these  modes. 

The  driving  force  behind  higher  frequency  and  therefore  higher  mode 
operation  was  the  feeling  that  high  frequency  operation  was  desirable  because 
better  collimated  beams  could  be  achieved,  thereby  reducing  the  effects  of 
adjacent  holes  In  the  resulting  scattered  fields.  However,  a  number  of 
disadvantages  readily  became  apparent:  (1)  the  operating  frequency  and 
excitation  efficiency  are  a  strong  function  of  plate  thickness  (member  B  In 
Fig.  1)  for  the  higher  order  modes  (n  >  0);  (2)  the  presence  of  bends  and  jogs 
affects  the  higher  order  modes  differently;  and  finally,  but  most  Importantly, 
(3)  the  actual  structure  was  found  to  exhibit  waveguldlng  properties  wherein 
energy  became  trapped  In  the  lower  lap,  and  propagated  along  the  Joint 
scattering  energy  from  many  holes.  This  made  It  Impossible  to  eliminate  the 
effects  of  adjacent  holes  with  the  frequency  limitations  Imposed  by  practical 
EMATs.  This  waveguldlng  effect  was  established  with  the  mock-up  which  became 
available  near  the  end  of  the  program  and  actually  proved  the  original 
transducer  configuration  to  be  unfeasible. 

Thus,  a  different  approach  became  necessary  and  the  original  far- 
fleld  scattering  Ideas  were  discarded.  First,  we  assembled  a  set  of 
transducers  having  an  0.5  In.  (1.27  cm)  period  which  could  be  operated  In  the 
n  «  0,  SH  mode  at  244  kHz,  thereby  making  the  operating  frequency  independent 
of  plate  thickness  and  the  excitation  efficiency  variation  a  weaker  function 
of  transducer  liftoff  and  plate  thickness.  Second,  we  located  the  transducer, 
as  Illustrated  in  Fig.  8,  as  close  as  was  feasible  to  the  hole  In  question, 
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thereby  performing  scattering  measurements  In  the  near  field.  Third,  we 
limited  the  transmitter  pulse  duration  to  a  2-5-cycle  tons  burst  and  gated  the 
received  signal  to  reject,  as  best  we  could,  the  direct  fsedthrough  signal  as 
well  as  signals  scattered  from  adjacent  holes  and  the  end  of  member  C 
(Fig.  8).  Fourth,  we  Implemented  a  measurement  system,  illustrated 
schematically  In  Fig.  9,  and  described  In  detail  In  section  v.  Fifth,  we 
viewed  the  section  A-A,  indicated  In  Fig.  8,  as  an  acoustic  bandpass  filter 
structure  having  a  predictable,  and  reproducible,  phase  and  amplitude 
response.  Thus,  by  taking  this  new  viewpoint,  the  inspection  scheme  becomes 
one  of  displaying  the  complex  Fourier  transform  of  the  received,  time-gated 
signals  which  Is  the  product  of  frequency  response  of  a  linear  system  composed 
of  the  structure,  EMAT  and  electronics.  The  presence  of  a  crack  modifies  the 
passband  shape  of  the  structure  and  thus  produces  a  change  In  the  Fourier 
transform.  We  are  thus  effectively  measuring  a  crack  Induced  change  In  the 
mechanical  response  of  the  joint. 

A  second  set  of  more  recent  experiments  on  actual  wing  sections  In 
which  the  fasteners  and  sealant  were  present  have  suggested  an  alternate, 
alblet  more  complicated,  measurement  scheme  Illustrated  In  Fig.  10.  Here 
three  transducers  are  placed  on  top  of  the  wing,  with  two  on  one  side  of  the 
front  and  a  third  on  the  other  side  so  that  the  transmission  properties  of  the 
n  =  0  SH,  mode  through  the  joint  can  be  taken  advantage  of.  Before  the  real 
wing  sections  became  available  we  reasoned  that  the  low  acoustic  Impedance  of 
the  sealant  wou  d  cause  the  two  parts  of  the  lap  joint  to  be  essentially 
uncoupled,  acoustically  so  that  this  transmission  would  not  occur.  This  Is, 
in  fact,  true  for  the  case  of  ultrasonic  waves  normally  Incident  upon  a  joint 
at  frequencies  of  several  MHz  as  would  be  the  case  if  a  piezoelectric 
transducer  were  placed  upon  the  top  surface  of  the  wing  joint  and  an  attempt 
were  made  to  make  a  high  frequency  inspection  of  the  lower  half  (see 
Appendix  A).  However,  for  the  n  =  0,  SH  waveguide  mode  propagating  along  the 
wlnq,  the  effect  of  the  sealant  is  negligible  at  low  frequencies,  as  discussed 
In  Appendix  A.  Thus,  for  the  geometry  of  Fig.  10,  an  n  *  0,  SH  wave  will 


16 


travel  across  the  joint  quite  efficiently.  This  was  a  completely  unexpected 
result,  hence  little  experimental  data  Is  available.  The  significance  of  the 
Fig.  10  configuration  is  that  independent  Information  Is  obtained  which  may 
complement  that  of  the  previous  configuration  to  help  discriminate  against 
effects  due  to  geometrical  changes. 


IV.  THEORETICAL  ANALYSIS  OF  SYSTEM 


The  configuration  shown  In  Fig.  8  was  selected  on  the  basis  of 
measurements  of  ultrasonic  scattering.  Fig.  11  compares  the  measured  SH 
wave  radiation  pattern  for  a  four-period  long,  n  =  0  EMAT  and  that  produced 
when  the  beam  strikes  an  0.25  In.  (0.635  cm)  dla.  hole.  The  data  was  obtained 

on  a  large  flat  plate  having  a  thickness  t  =  0.25  in.  (0,635  cm),  by  locating 

the  transmitter  the  same  distance  from  the  hole  as  In  Fig.  8  and  rotating  the 
receiver  about  the  hole  while  keeping  it  aimed  at  the  hole.  For  comparison, 
we  have  illustrated  the  wing  section  member  E  of  Fig.  8  by  dashed  lines  at  the 
location  that  it  would  occupy  In  the  proposed  Inspection  configuration. 

It  should  be  noted  that,  had  we  rotated  the  transducer  about  the 
center  of  the  transmitter,  a  somewhat  narrower  beam  pattern  would  have  been 
measured  in  which  the  -3  db  width  was  46°  and  the  -6  db  width  was  60°. 

The  key  results  deduced  from  these  radiation  pattern  and  scattering 
measurements  are:  (1)  energy  is  confined  about  a  single  hole  In  the 

Inspection  scheme  illustrated  in  Fig.  8  and  (2)  on  forward  transmission  the 

hole  modifies  the  radiation  pattern  of  the  transducer,  but  considerable  energy 
(only  .9  db  forward  loss)  passes  by  the  hole  to  reflect  from  the  edge  of 

member  E  in  Fig.  8,  and  (3)  little  energy  is  scattered  at  90°.  Thus,  when  the 

transducers  are  positioned  as  shown,  the  hole  can  be  expected  to  have  minimal 
influence  on  the  ultrasonic  signals.  Based  upon  these  radiation  pattern 
measurements,  we  can  derive  a  simple  model  of  the  ultrasonic  reflection 
process  in  the  absence  of  sealant  and  fasteners;  i.e.,  only  the  part  of  the 
lap  joint  shown  in  Fig.  8. 

Fig.  12(a)  is  a  schematic  representation  of  the  entire  signal 
transmission  process.  A  two-cycle,  244  kHz,  tone  burst  is  input  to  a  four- 
period  EMAT,  thereby  exciting  an  n  =  0,SH  wave  having  an  0.5  in.  (1.27  cm) 

long  wavelength.  The  SH  wave  is  filtered  by  the  wing  section  geometry  of 


20 


SIMPLE  WING  SECTION  in  l 

EQUIVALENT  CIRCUIT 

<b)  SIMPLIFIED  LINEAR  EQUIVALENT  CIRCUIT  DUE  TO  TIME  GATING. 


Fig.  8  and  subsequently  received  by  an  identical  receiving  EMAT,  amplified  and 
filtered.  As  indicated  in  Fig.  12(a),  each  element  in  the  chain  has  a 

frequency  response. 

The  objective  is  to  measure  the  effect  a  crack  has  on  the  wing 
section  transfer  function  Fs((u),  and  thereby  deduce  the  length  of  the  crack. 
The  simplified  equivalent  circuit  of  Kig.  12(b)  is  arrived  at  as  follows. 

HU)  denotes  an  equivalent  system  impulse  response  which  is  determined  by 
calculation  or  by  measuring  the  direct  transmission  of  an  SH  wave  between  two 
transducers  on  a  large  9/32  in.  (0.71  cm)  plate.  Time  gating  is  used  to 
eliminate  the  eff.cts  of  electromagnetic  feedthrough  and  other  Interfering 
signals  such  as  reflections  from  sample  edges.  In  addition  to  the  EMAT 
response,  FjU),  the  spectral  characteristics  of  the  driving  voltage,  V ^ n ( uj ) , 
and  of  the  receiving  amplifier,  RU),  are  included. 

The  actual  wing  section  transfer  function  Fs(«)  is  extremely 
difficult  to  calculate  since  it  is  due  to  a  very  complicated  and  analytically 
intractable  structure.  However,  if  we  restrict  the  frequency  band  of  Interest 
by  designing  HU)  appropriately,  then  based  on  the  scattering  data  discussed 
above  as  well  as  experimental  observation  a  first  order  approximation  is  to 
replace  FSU)  by  S0U)  as  in  Fig.  12(b).  Here  we  assume  that  energy  can 
propagate  from  transmitter  to  receiver  by  either  oF  the  two  paths  shown  in 
Fig.  8.  Path  1  is  due  to  a  reflection  off  the  stpp  of  Fig.  8,  characterized 
by  an  amplitude  transfer  function  yU)  and  a  delay  path  21^.  The  second 
signal  path  has  a  length  2L2  and  an  amplitude  response  is  characterized  by 
d(oj) •  This  corresponds  to  energy  entering  the  joint  and  reflecting  from  its 
edge.  aU)  Is  what  we  are  seeking  because  it  is  affected  by  the  presence  and 
shape  of  the  fastener  hole. 

The  reason  there  Is  not  a  third  delay  path  corresponding  to  the  hole 
location  is  due  to  the  fact  that  little  energy  is  scattered  at  90°  from  a  hole 
as  was  illustrated  in  Fig.  11. 
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An  actual  wing  section  win  exhibit  variations  In  the  thickness  of 
members  B  and  C,  Fig.  8,  which  affects  both  HU)  and  SQU),  as  well  as 
''■aria' ‘oiij  in  the  length  ui  member  C  which  affects  and  aU).  Furthermore, 
the  presence  of  the  sealant  win  cause  energy  to  be  coupled  to  the  second 
member.  Even  so  there  will  still  be  two  main  signal  paths,  and  L£.  due  to 
the  discontinuities  indicated  In  Fig.  12.  There  will  be  additional 
information  about  the  possibility  of  flaws  in  the  upper  half  of  the  joint 
because  energy  Is  coupled  through  the  sealant  Into  this  region  as  well. 

To  effectively  deal  with  these  structural  variations,  a  calibration 
procedure  must  be  arrived  at.  With  the  aid  of  Fig.  12(b)  we  see  that  l_j  Is 
set  by  the  transducer  placement  and  Is,  therefore,  known.  Lg  will  be  known 
from  the  wing  design  and  yU)  will  change  in  a  simple  and  controlled  way 
because  the  step  height  does  not  change  that  much.  The  success  of  the 
Inspection  process  will  hinge  upon  associating  changes  in  the  transmitted 
slgn.il  1  i  changes  in  dU),  Independent  of  these  known  varatlons. 

In  order  to  develop  a  model,  we  will  calculate  HU)  and  SQU) 
independently.  Vaslle  and  Thompson^’2'  give  expressions  for  SH  wave 
excitation  in  metal  plates.  For  present  purposes,  It  is  possible  to  simplify 
and  adapt  these  results  In  order  to  exhibit  some  Important  features  that  have 
been  observed  experimentally.  Thus,  by  assuming  zero  skin  depth,  neglecting 
transducer  end  effects  and  diffraction,  and  noting  that  the  forward  transfer 
Impedance  1$  proportional  to  u/Bn  Un  ■  SH  plate  mode  propagation  constant)  we 
can  derive  an  expression  for  the  frequency  dependence  of  HU): 

HU)  »  Hfl.  RU)  •  TU)  •  |zoU)  +  Z^m)  +  etc.]  (1) 

where  H0  is  an  arbitrary  scale  factor,  RU)  is  the  electronic  receiver 
frequency  response,  T(u>)  is  the  Fourier  transform  of  the  Input  current 
waveform  given  approximately  by: 

TU)  -  sine  [i, ( f-fQ )  Nc/fQ]  (2) 
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a)  =  2itf,  sine  x  =  sin  x/x,  Nc  =  number  of  carrier  cycles,  f0  =  carrier 
frequency  and  the  bracketed  terms  give  the  transducer  response  denoted  F^w) 
In  Fig.  12.  This  analysis  differs  slightly  from  that  shown  In  Fig.  12(b) 
since  we  are  specifying  the  Input  current,  Km),  rather  than  voltage  V^n(w). 


The  bracketed  terms  In  Eq.  (1)  are  expressed  as  a  sum  of  normalized 
forward  transfer  Impedances  of  the  plate  modes  which  are  dependent  upon  the 
EMAT  parameters,  plate  material  and  thickness.  For  present  purposes  we 
consider  the  first  two  modes  and  have: 


Similar  terms  would  be  obtained  for  ZnU),  n  >  1.  Here:  V$  =  shear  wave 
velocity  (3111  m/s),  LD  Is  the  length  of  the  transmission  path,  f0  =  (Vs  * 
transducer  period),  fcl  ■  (nVs/2t)  nsl  =  218  kHz  =  cutoff  frequency  of  n  »  1 
mode  for  the  t  =  9/32  in.  plate  thickness,  t  =  9/32  In.,  and  N  a  number  of 
periods  In  the  transducers  =  4. 

The  reason  for  Including  only  the  first  two  modes  Is  demonstrated  by 
the  dispersion  relation  shown  in  Fig.  13.  The  long  dashed  line  corresponds  to 
a  set  of  0.5  In.  (1.27  cm)  period  EMATs  operating  on  a  9/32  In.  (0.71  cm) 
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Fig.  13  SH-Dlsperslon  Relation  for  a  9/32"  Aluminum  Plate. 


plate.  At  any  point  where  the  dashed  line  crosses  a  solid  line,  a  single  inode 
can  be  excited  by  a  monochromatic  drive.  For  a  tone  burst,  several  modes  can 
be  excited  with  amplitudes  determined  by  the  spectral  content  of  the  burst. 

For  our  case, the  n  =  0  and  n  ■  1  modes  are  important. 

The  transducers  ultimately  used  were  not  originally  designed  for  this 
project  and  became  available  late  in  the  program  from  other  work.  Hence  they 
were  not  fully  optimized  for  this  particular  problem.  Because  of  time 
resolution  problems  associated  with  electromagnetic  feedthrough  suppression 
and  geometrical  reflections,  we  limited  the  lengths  of  the  EMATs  to  four 
periods  (2  In.  (5  cm)  long)  and  limited  the  250  kHz  drive  current  burst  to  two 
or  three  cycles,  even  though  eight  cycles  would  have  provided  maximum  signal 
levels  as  well  as  better  higher  order  mode  suppression.  The  consequence  of 
these  limitations  are  Illustrated  In  Fig.  14,  which  contains  Individual 
sketches  of  RU),  TU),  Zq{ u> )  ,  and  ZjU)  for  fc  «  250,  N  ■  4,  and 
Nc  ■  2  In.  (5  cm).  Also  shown  are  sketches  of  the  combination  appearing  In 
Eq.  (1)  which  seem  to  give  HU).  From  the  lower  plot  It  Is  clear  that  the 
n  ■  1  mode  will  contribute  significantly  In  a  direct  transmission  measurement 
and  this  has  been  observed  experimentally,  Thus  a  better  system  design 
requires  lower  frequency  operation  and  longer  period  transducers  so  that  the 
n  »  1  mode  Is  cutoff  and  Z^U)  <<  ZqU)  over  the  band  of  Interest. 

Fortunately,  for  the  geometry  of  Fig.  8  where  the  acoustic 
transmission  path  Involves  a  complicated  reflection/scattering  process,  It 
appears  that  the  contribution  of  ZjU)  Is  reduced  because:  (1)  the  7/32  In. 
(0.56  cm)  thick  member  E  Is  cutoff  to  the  n  *  1  mode  at  frequencies  below 

280  kHz  so  that  It  does  not  contribute  to  path  2  and  (2)  the  Incident  n  *  1 

mode  appears  to  Incur  a  higher  loss  upon  reflection  off  the  step  so  that  Its 

effect  Is  lower  In  path  1  than  It  is  In  direct  transmission.  Thus,  based  upon 

these  observations,  It  Is  reasonable  to  neglect  the  n  ■  1  mode  In  developing  a 
first  order  model  of  the  wing  structure  for  the  n  =  0  mode  as  was  done  In 
drawing  the  equivalent  circuit  In  Fig.  12(b). 
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Fig.  14  Components  for  H  (u)  for  Direct  Transmission. 
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From  the  single  mode  equivalent  circuit  of  Fig.  12(b)  we  compute  the 
structure  transfer  function  SqU)  which,  when  multiplied  by  an  appropriate 
transfer  function  HU),  gives  the  overall  system/transfer  function 
S0U)  x  HU)  which  we  measure.  The  result  Is: 


or 
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From  measurements  on  the  geometry  of  the  wing  inock-up  Illustrated  In 
Fig.  8:  =  1.85  In.  (4.7  cm);  2L^/VS  =  30.33  us;  L.?  *  2.65  In.  (6.7  cm); 

2I_2  »  43.44  ys;  At  ■  2L2/VS  -  2L j /V$  3  13.11  ps.  Measurements  suggest  that 
yU)  Is  a  slowly  varying  function  of  frequency  as  Is  aU).  Thus,  It  Is  useful 
to  examine  the  behavior  of  (6)  with  5U)/yU)  3  A,  a  constant.  Thus: 


|  SqU)/2y  (u>)  |  »  cos  ( tt f at  ) ;  A  3  1 


(7) 


3  1.25  +  cos  (2tt fAr) ;  A  ■  .5 


(8) 


Fig.  15(a)  Is  a  plot  of  Eq.  (7)  and  (8),  from  D.C.  to  400  kHz. 

Because  the  n  »  1  mode  Is  cutoff  In  member  E  of  Fig.  8,  a  periodic 
response  similar  to  that  of  Fig.  15  does  not  exist. 

Fig.  15(b)  combines  the  results  presented  in  Fig.  14  with  15(a), 
neglecting  the  Z^U)  In  HU),  and  renormal izlng  the  results  to  unity  at 
250  kHz. 
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(cj)  =  R(w)-T|w)  *  }ZJ(J)| 


The  Importance  of  Fig.  15(a)  is  that,  for  A  =  1,  one  would  expect  to 
observe  a  characteristic  null  in  the  Fourier  transform  of  the  system  Impulse 
response  near  265  kHz  corresponding  to  destructive  cancellation  of  SH  signals 
traveling  along  the  two  delay  path  illustrated  In  Fig.  12.  If  an  obstruction 
is  Introduced  into  the  L2  path,  e.g.,  a  crack  emanating  from  the  fastener 
hole,  then  the  null  should  be  strongly  affected  as  indicated  in  Fig.  15(b)  for 
A  *  .5.  Both  theoretical  predictions  have  In  fact  been  observed 
experimentally  and  will  be  described  below.  Before  concluding  the  section  It 
should  be  pointed  out  that  the  structure  response  exhibits  characteristic 
features  that  land  on  the  skirts  of  our  EMAT  test  system  response  where  the 
signal  to  noise  ratio  Is  reduced.  Improved  performance  could  be  achieved  by 
shifting  the  EMAT  passband  down  to  196  kHz.  These  numbers  would,  of  course, 
be  different  for  samples  with  dimensions  different  from  those  of  the  mock-up 
used  in  this  work. 

Early  effort  was  directed  towards  using  .25  In.  (.64  cm)  period  EMATs 
at  higher  frequencies  (650-850  kHz).  Referring  to  the  SH  mode  dispersion 
relation  in  Fig.  13,  for  a  9/32  In.  (0.71  cm)  plate  t/x  «  1.125  and  operation 
with  a  single  higher  order  SH  mode  is  possible  for  n  ■  2  or  n  ■  3  by 
restricting  the  system  bandwidth.  Such  an  approach  offered  several  advantages 
such  as  Improved  signal  to  noise  ratio,  better  time  resolution  due  to  the 
slower  velocities,  and  better  beam  collimation  because  of  high  frequency 
operation.  The  major  disadvantage  was  that  frequency  of  operation,  excitation 
efficiency,  and  group  delay  were  strong  functions  of  plate  thickness.  Since 
actual  wing  sections  do  have  significant  geometrical  variations,  such  added 
variability  would  further  complicate  data  Interpretation.  Also  the  sealant  In 
actual  wing  sections  would  affect  the  higher  order  SH  modes  In  a  more 
uncontrolled  manner. 


V.  DESCRIPTION  OF  MEASUREMENT  SYSTEM 


A.  ELECTRONICS 

A  series  of  breadboard  electronics  were  partially  developed  and 
Improved  on  the  program  In  order  to  facilitate  data  Interpretation  and  pattern 
recognition  In  real  time.  Figure  9  gives  a  block  diagram  of  the  measurement 
system  comprising: 

Function  Generator  -  A  breadboard  function  generator  supplies  digital 
drive  signals  to  the  a  class  C  power  amp.lfler  at  a  low  pulse 
repltltlon  rate.  The  drive  signal  period,  prf  and  pulse  durations 
are  all  adjustable.  Fig.  16  shows  the  transmitter  drive  voltage, 
resulting  EMAT  current  waveforms  as  well  as  the  received  direct 
transmission  acoustic  signals  for  two  cases.  There  Is  significant 
distortion  present  In  the  current  because  the  transmitter  was 
designed  for  operation  at  750  kHz  and  simply  modified  for  250  kHz 
operation.  In  the  first  case,  an  8-cycle,  250  kHz,  tone  burst  Is  the 
Input,  corresponding  to  a  maximum  amplitude  received  signal.  (The 
current  rundown  Is  by  design  to  protect  the  transmitter  from  long 
duration  pulses,  since  no  heat  sinking  is  used.)  The  second  shows 
the  effect  of  reducing  the  drive  to  two  TTL  pulses  as  Is  necessary  In 
order  to  minimize  Interference  from  the  direct  electromagnetic 
feedthrough  pulse.  As  the  EMATS  act  as  bandpass  filters,  the  current 
waveform  distortion  Is  of  little  practical  significance.  However, 
the  fact  that  a  strong  second,  higher  frequency  trailing  pulse 
appears  when  the  Input  pulse  duration  Is  reduced  to  two  carrier 
cycles  clearly  shows  the  existence  of  the  n  ■  1  mode  which  propagates 
slower  and  contains  more  high  frequency  energy  than  the  n  ■  0  mode. 

The  pulse  repltltlon  rate  used  was  40  pps,  or  25  ms/pulse,  a 
compromise  figure.  Higher  rep-rates  coupled  with  low  pass  filtering 
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M  I)  SKiNAl. 

Fit],  16  Ultrasonic  Waveforms. 

Top:  TTL  Transmitter  Drive  and  EMAT  Current  Waveforms 

(5  Amp/ Di v )  for  8  Cycle  (Left)  and  2  Cycle  (Right) 
Exc i tat i on . 

bottom:  Received  Ultrasonic  Signal  on  a  9/32"  Aluminum 

Plate.  The  Richer  Spectral  Content  of  the  Shorter 
Pulse  Excites  the  n  =•  1  Mode. 
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is  advantageous  for  averaging  out  thermal  noise  and  rapid  data 
acquisition.  However,  since  the  acoustic  Q  of  the  aluminum  mock-up 
was  very  high  at  250  kHz,  reverberations  did  not  decay  rapidly.  Thus 
a  high  rep-rate  coupled  with  time  jitter  led  to  significant  noise 
like  clutter  interference.  Fortunately,  the  actual  wing  section  was 
long  and  narrow  and  appeared  to  have  a  significantly  lower  acoustic  Q 
(possibly  due  to  the  sealant).  Hence  this  type  clutter  interference 
was  not  a  problem. 

Receiver  -  The  receiver  shown  in  Fig.  9  Is  a  four-stage  discrete, 
fa^t  recovery  FET  amplifier  capable  of  80  dB  voltage  gain.  Simple 
single  tuned  stages  give  the  response  Rv •  ■ )  illustrated  in  Fig.  14. 


Correlation  Processor  •  The  analog  correlation  processor  forms  the 
linear  product  of  the  received  acoustic  signals  with  a  known,  gated 
reference  signal  in  order  to  derive  both  phase  and  amplitude 
Information  about  received  acoustic  signals.  By  using  the  VCO 
capability  of  the  HP3312B,  it  is  possible  to  sweep  the  frequency  of 
the  gated  reference  and  produce  a  real  time  plot  corresponding  to  the 
Fourier  transform  of  tine  gate  segments  of  the  received  acoustic 
signals.  Fig.  17  gives  functional  block  diagram  of  the  processor 
board.  The  Ref.-l  and  Ref. -2  Inputs  are  identical  TTL  signals,  half¬ 
bit  delayed  with  respect  to  each  other,  and  time  delayed  with  respect 
to  the  main  transmit  TTL  burst.  The  reference  signals  are  buffered 
and  bandpass  filtered  in  order  to  suppress  higher  harmonics  which 
degrade  signal  to  noise.  Thus, the  ideal  output  of  the  processor  may 
be  approximated  by: 


V^R1  = 


v  (t) 

a 


C  OStiin 
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Processor  Functional  Block  Diagram. 


W  =  C<I  2  va(t)  s1nURtdt> 

'l 

where  C  Is  a  constant  of  proportionality  and  is  controlled  by  the 
receiver  gain  adjust.  The  brackets  denote  time  average.  If  va(t)  Is 
a  signal  that  vanishes  outside  the  range  T ^  <  t  <  T2  then  Vx  and  Vy 
are  the  real  and  Imaginary  parts  of  Fourier  transform,  /va(t), 
evaluated  at  w  =  up.  Since  the  complex  Fourier  transform  has  fewer 
wiggles  than  the  corresponding  time  waveform  va(t)  In  the  case  of 
overlapping  time  signals,  It  Is  a  much  simpler  display  to 
Interpret.  The  time  averaging  Is  achieved  by  low  pass  filtering  the 
sample  and  hold  outputs.  As  a  comprorlse  between  speed  of  data 
acquisition  and  degree  of  averaging  we  used  a  single,  two-stage  RC 
low  pass  filter  and  set  each  time  constant  at  .033  sec.  Thus,  we 
were  able  to  effectively  average  about  seven  sample  and  hold  outputs 
at  a  40  pps  rep-rate  (e.g.,  Neff;  *  (2  x  .033  x  40) 2 ) .  The  price  paid 
for  averaging  Is  a  slower  data  acquisition  rate. 

0.  DATA  DISPLAY 

Two  modes  of  display  were  used.  Polar  plots  of  the  locus  of  the 
complex  Fourier  transform  were  obtained  by  using  the  Vx  and  Vy  outputs  to 
drive  an  x-y  recorder.  Plots  of  the  magnitude  of  the  filter  passband  were 
obtained  from  the  VCO  voltage  and  log  outputs.  Variation  of  both  the 
transducer  orientation  (angular  dependence  of  scattering)  and  reference 
frequence  (Fourier  transform)  were  considered.  Because  of  geometrical 
variation  In  the  actual  structure,  more  repeatable  data  was  obtained  from  the 
Fourier  transform  approach  with  fixed  transducer  location  and  variable 
frequency. 
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The  system  dynamic  range  was  limited  to  about  25  dB  by  error  in  the 
squaring  circuit.  Thermal  noise  was  found  to  be  below  the  clutter  level  due 
to  acoustic  reverbation  in  the  wing  section  mock-up.  Excellent  signal  to 
noise  was  observed  in  the  actual  wing  section  because  its  long  and  narrow 
shape  as  well  as  damping  Introduced  by  the  sealant-filled  lap  joint  and 
fasteners  reduced  the  clutter. 

Fig.  18  shows  the  various  types  of  displays  for  a  direct  trans¬ 
mission  SH  signal  on  a  9/32  In.  (0.71  cm)  thick  plate.  The  top  photograph 
shows  the  received  time  signal  In  relation  to  the  "2  cycle"  TTL  transmitter 
drive  burst  superimposed  on  one  of  the  TTL  reference  signals.  The  power 
spectra  (c),  plotted  from  200-300  kHz  clearly  shows  the  presence  of  the  n  -  1 
mode  as  predicted  above  and  illustrated  in  Fig.  14.  Part  (d)  is  a  polar 
plot  of  the  complex  Fourier  Transform.  The  fluctuation  observed  in  (c)  and 
(d)  are  not  due  to  thermal  noise,  but  to  acoustic  reverberation  as  discussed 
above. 


C.  EXPERIMENT  rETUP 

Fig.  19  .  .  Ketch  of  the  wing  section  mock-up  we  fabricated  in¬ 
cluding  the  transducer  locations.  This  sample  was  produced  in  our  machine 
shop  to  serve  as  a  basis  for  experimentation  until  real  samples  could  be 
obtained.  Fig.  20  Is  a  photograph  of  the  experimental  setup  showing  the 
transducer  fixture  In  position  on  a  wing  section  supplied  later  by  the  Air 
Force.  The  transducers  are  held  in  fixed  positions  with  respect  to  each  other 
by  the  aluminum  jig.  The  pin  is  used  as  a  marker  to  position  the  apparatus 
properly  with  respect  to  the  fastener.  The  breadboard  custom  electronics  is 
seen  toward  the  rear  of  the  setup. 
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Fig.  18  Time  Waveforms  and  Real  Time  Fourier  Transforms. 


VI.  EXPERIMENTAL  RESULTS 


Four  sets  of  experiments  were  performed.  During  the  early  part  of 
the  program,  attention  was  directed  toward  the  higher  order  modes.  These 
measurements  were  made  on  the  wing  suction  mock-up  sketched  In  Fig.  19,  As 
discussed  previously,  their  sensitivity  to  geometrical  changes  In  many 
features  of  the  structure  made  tht  data  obtained  unreliable,  and  hence  results 
are  not  reproduced  here.  As  a  result,  the  fundamental  mode  was  Investigated 
on  the  mock-up,  as  discussed  In  section  A  below.  Analysis  of  this  data  led  to 
the  theoretical  model  presented  In  section  IV.  Toward  the  end  of  the  program, 
a  section  of  a  real  wing  became  available.  Section  B  presents  the  preliminary 
data  obtained  on  this  sample. 

A.  MOCK-UP 

Figs.  8  and  9  Illustrated  the  experimental  setup  in  which  only  the 
lower  section  of  the  lap  joint  was  tested  In  the  absence  of  fasteners  and 
sealants.  The  sample  was  first  used  to  optimize  the  electronics  and  establish 
a  measurement  technique  that  assured  repeatability  for  an  unperturbed  set  of 
holes.  Then  slots  were  placed  In  several  holes  with  a  jeweler's  saw.  The 
slots  were  cut  toward  the  outer  edge  of  the  lower  lap  as  was  Illustrated  In 
Fig.  19. 


The  experimental  results  obtained  on  the  mock-up  are  given  In 
Fig.  21.  While  the  time  waveforms  are  difficult  to  Interpret,  the  power 
spectra  are  not.  The  sharp  dip  appearing  In  the  power  spectrum  at  270  KHz 
for  an  unslotted  hole,  Fig.  21(a),  Is  In  gona  agreement  with  tho  simple  model 
presented  above.  The  frequency  of  the  dip  Is  at  the  predicted  location  and, 
as  the  slot  length  Is  Increased,  the  amplitude  of  the  dip  Is  decreased.  This 
closely  corresponds  to  the  behavior  predicted  in  Fig.  15.  Also  shown  In 
Fig.  21  are  the  complex  Fourier  transform  plots.  From  these  It  can  be  seen 
that  there  Is  phase  as  well  as  amplitude  Information  that  varies 
systematically  with  slot  length.  These  correlations  have  not  yet  been  fully 
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Progression  on  Wing  Section  Mockup 


pursued.  Fig.  22  demonstrates  the  correlation  between  slot  length  and 
experimental  results  as  derived  from  only  the  power  spectra  plotted  in 
Fig.  21.  These  results  were  quite  repeatable  and  encouraging.  In  par¬ 
ticular,  the  depth  of  the  dip  at  270  kHz  varies  monotonlcally  with  slot 
length.  The  other  features  provide  information  about  the  geometry  that  has 
not  yet  been  fully  explored. 

B.  WING  SECTION 

After  establishing  that  a  reliable  technique  had  been  arrived  at  In  a 
controlled  experiment  on  a  geometrically  uniform  wing  section  mock-up  with  no 
sealant  or  upper  half  of  the  joint,  we  proceeded  to  evaluate  the  technique  on 
an  actual  wing  section  In  which  both  the  fasteners  and  sealant  were  present. 

No  cracks  were  known  to  be  present  In  any  of  the  holes  at  the  beginning  of  the 
experiment.  The  geometry  of  this  wing  section,  shown  previously  In  the 
photograph  of  Fig.  20,  Is  Illustrated  In  Fig.  23,  It  can  be  seen  that  there 
are  regions  in  which  both  the  width  of  the  lap  joint  and  the  thickness  of  the 
plate  vary.  All  holes  are  not  shown,  but  those  at  points  where  the  geometry 
changes  are  included. 

Fig.  24  shows  a  set  of  data  obtained  for  34  fasteners  in  this  wing 
section.  The  data  Is  cleaner  than  that  obtained  onthe  mock-up  and  presented 
In  Fig.  21  because  of  the  lower  acoustic  Q  of  this  sample.  Again  a  dip  Is 
observed.  However,  It  Is  shifted  down  In  frequency  because  of  an  Increase  In 
lap  width  with  respect  to  that  of  the  mock-up.  The  frequency  of  240  kHz  Is 
again  In  agreement  with  the  predictions  of  the  model. 

Inspection  of  Fig.  24  shows  that  there  are  considerable  differences 
In  the  responses  of  the  Individual  holes,  and  this  can  be  related  to  the 
changes  In  sample  geometry  shown  In  Fig.  23.  Note,  however,  that  curves  15 
through  25  exhibit  extremely  uniform  responses  because  the  wing  section  Is 
uniform  In  this  region.  This  data  then  demonstrates  the  fact  that 
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Feature  Variation  with  Crack  Depth 


Geometry  of  Lap  Joint 


reproducible  data  can  be  obtained  on  real  wing  sections,  but  that  geometrical 
changes  also  influence  the  data. 

These  measurements  show  the  background,  "good  fastener,"  response, 
which  Is  somewhat  different  from  that  shown  in  Fig.  21  for  the  mock-up  because 
the  increase  in  width  of  the  lap  shifted  the  Interface  null  down  In  frequency, 
and  the  presence  of  the  sealant  and  upper  half  of  the  joint  changes  the 
parameters  model ed  by  5  and  y  and  thus  changes  the  envelope  of  the  filter 
passband.  To  test  if  the  same  procedure  for  detecting  cracks  would  work  in 
this  case,  we  performed  a  controlled  experiment  wherein  we  removed  a  fastener 
and  then  proceeded  to  cut  an  0.100  in.  (.254  cm)  slot.  We  were  able  to  detect 
the  presence  of  this  slot,  as  summarized  in  Fig.  25.  Two  points  are 
illustrated  by  the  data.  First,  the  presence  or  abienre  of  the  bolt  does,  in 
fact,  change  the  experimentally  observed  spectrum.  Second,  under  constant 
geometrical  conditions,  with  bolt  removed,  the  presence  of  a  crack  in  the  he1'- 
affects  the  spectrum  In  the  same  way  as  before.  By  comparing  the  last  two 
sets  of  data,  one  sees  that  the  depth  of  the  null  decreases  as  a  crack 
appears,  just  as  was  shown  for  the  mock-up  in  Fig.  22. 

A  final  set  of  experiments  was  performed  in  the  transmission  geometry 
illustrated  in  Fig.  10.  Here,  measurements  were  made  of  the  amplitude  of  the 
ultrasonic  signal  transmitted  from  one  side  of  the  lap  joint  to  the  other.  It 
is  Initially  surprising  that  this  transmission  occurs  In  light  of  previous 
experiments  In  which  ultrasonic  energy,  injected  through  the  top  half  of  the 
lap  joint,  did  not  effectively  penetrate  through  the  sealant  into  the  lower 
half.  Reasons  for  this  difference  are  discussed  in  Appendix  A. 

Fig.  26  shows  a  waveform  photo  and  the  power  spectra  data  for 
the  geometry  of  Fig.  10  in  which  the  transmission  through  the  joint  was 
monitored.  Thus  the  wave  was  launched  by  the  upper  left  transducer  in 
Fig.  10,  and  was  detected  by  the  lower  right  transducer  after  passing  through 
the  joint. 
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Transmission  Power  Spectra 


A  rather  dramatic  change  in  the  energy  transmitted  past  the 
previously  slotted  hole  No.  19  and  Into  the  other  half  of  the  lap  joint  is 
observed,  even  after  the  fastener  had  been  reinstalled.  It  appears  likely 
that  the  presence  of  the  crack  substantially  Interrupts  the  transmission  of 
energy  across  the  Interface  and  thereby  Is  being  accurately  detected,  a  quite 
exciting  result.  However,  further  experiments  are  needed  to  be  certain  that 
the  result  was  not  caused  by  damage  to  the  sealant  during  the  removal  and 
reinsertion  of  the  fastener. 


VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  five-month  feasibility  study  we  have  demonstrated  that: 

(1)  EMATS  placed  on  the  outer  wing  surface  can  be  used  to  excite  ultrasonic 
energy  which  can  propagate  to  and  return  from  the  fastener  region  of  the  lower 
section  of  a  wing  lap  joint,  (2)  Information  relating  to  the  fastener  hole 
integrity  is  contained  In  both  signals  reflected  from  the  joint  and 
transmitted  through  the  joint  to  the  other  plates,  and  (3)  real  time  Fourier 
analysis  Is  useful  in  deducing  flaw  Information  from  the  complex  ultrasonic 
response.  For  the  case  of  a  lap-joint  of  controlled  geometry  and  without 
sealant  or  fastener  present,  It  has  been  established  that:  (4)  saw  slots  as 
small  as  0.030  in.  can  be  detected  and  (5)  an  acoustic  model  has  been 
developed  which  explains  many  features  of  the  experimental  data  and  aids  In 
its  1 nterpretatlon  and  In  system  optimization.  A  limited  numer  of  experiments 
has  been  performed  on  complete  lap  joints.  Encouraging  results  have  been 
obtained,  Including:  (6)  the  detection  of  a  0.10D  In  (.254  cm)  slot  and 
(7)  the  establishment  that  the  Instrumental  response  Is  reproducible  In 
regions  of  constant  geometry.  However,  further  work  Is  required  to  separate 
the  flaw  Information  from  that  Introduced  by  variations  In  overlap  of  the 
plates,  oy  the  fastener,  and  by  coupling  to  the  other  plate  through  the 
sealant  In  regions  where  these  geometrical  features  are  not  constant. 

Future  efforts  will  require:  (1)  optimization  of  system  design  to 
Improve  the  precision  of  the  data  (2 )  development  of  more  sophisticated  data 
gathering,  processing,  and  Interpretation  procedures  to  separate  the  flaw 
information  from  the  unwanted  background,  and  (3)  a  more  extensive,  low 
frequency  experimental  verification  to  establish  a  statistically  significant 
foundation  for  field  application  of  this  technique.  A  more  detailed  plan  will 
be  submitted  under  a  separate  cover. 
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APPENDIX  A  —  EFFECTS  OF  SEALANT 


Fig.  A-l  illustrates  two  cases  of  practical  interest  in  inspecting 
lap  joints:  (1)  that  of  a  shear  or  longitudinal  wave  normally  Incident  upon 
an  alumlnum-polyethelene-alumlnum  interface  and  (2)  that  of  an  n  s  0,  SH  mode 
propagating  along  an  infinite  aluminum-polyethelene-aluminum  sandwich 
structure.  Polyethelene  is  used  to  simulate  the  acoustic  effects  of  a 
sealant,  whose  properties  we  have  not  characterized.  Case  (1)  relates  to  the 
situation  which  occurs  during  Inspection  through  the  upper  sheet  with  a  device 
like  the  Rotoscanner.  Case  (2)  occurs  during  inspection  of  the  type  discussed 
in  this  report. 

Normal  Incidence 

The  respective  shear  and  longitudinal  wave  velocities,  in  units  of 
km/sec,  for  aluminum  (density  pa  =  2690  kg/m3)  and  polyethelene 
(pp  =  886  kg/m3)  are:  Vs  =  3040,  540;  VL  *  6420,  1950.  The  relative  acoustic 
Impedances  of  the  two  materials  are  Zp/Za  ■  .0585  and  .1011  for  shear  and 
longitudinal  waves,  respectively.  We  measured  the  thickness  of  the  sealant  to 
be  about  5  mils  in  our  sample.  It  possibly  varies  from  1-10  mils  throughout 
the  structure.  Taking  a  5  mil  nominal  thickness  we  find  that  the  adhesive 
thickness  is  X/4  at  1.06  MHz  for  a  normally  incident  shear  wave  and  at 
3.84  MHz  for  a  longitudinal  wave.  Thus  It  follows  from  simple  transmission 
line  analysis  that  strong  reflection  loss  peaks  will  occur  at  the  quarter  wave 
frequencies  given  by  (2  m  +  1)  x  1.06  MHz  and  (2  m  +  1)  x  3.84  MHz  for  shear 
and  longitudinal  waves,  respectively.  Thus,  for  a  double  pass,  a  loss  due  to 
reflection  of  36  dB  for  shear  waves  and  28  dB  for  longitudinal  waves  Is  to  be 
expected  at  their  respective  quarter  wave  frequencies. 

At  first  glance  one  would  expect  transmission  bands  to  occur  at  the 
half  wave  frequencies,  but  as  ■indicated  In  Fig.  A-2,  they  are  quite  narrow. 
Furthermore,  a  variation  in  sealant  thickness  of  from  5. 7-8. 8  mil  would  shift 
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Transmission  Through  0.005"  Thick  Polyethelene  Layer 


the  quarter  wave  frequency  of  a  shear  wave  transmission  over  the  3-6  MHz 
band.  One  would  then  expect  severe  phase  distortion  of  a  pulse  and  strong 
variation  from  point  to  point  on  the  sample.  Similar  effects  are  expected  for 
longitudinal  wave  transmission.  Hence,  It  Is  not  practical  to  Inspect 
ultrasonlcally  the  lower  lap  joint  from  the  top  surface  with  a  small,  high- 
frequency  piezoelectric  transducer  probe. 

Plate  Modes 

Experiments  with  low  order  SH  waves  In  sandwich  structures  of  the 
type  In  Fig.  A-l(b)  Indicate  that,  at  a  frequency  of  250  kHz,  the  n  *  0  mode 
Is  only  slightly  perturbed  by  the  sealant.  At  this  frequency,  the  thin,  slow 
layer  can  transmit  a  shear  stress  easily  from  the  top  plate  to  the  bottom 
plate.  Thus  good,  n  =  0,  SH  wave  transmission  Is  to  be  expected  through  the 
lap  joint,  with  reflection  Introduced  by  the  step  discontinuities.  Howev-r, 
sealant  variations  may  Introduce  uncontrolled  amplitude  variations  In  the 
transmission  through  and  reflection  from  a  lap  joint  If  the  frequency  Is  not 
sufficiently  low  enough  to  cut  off  all  higher  order  modes  that  can  be 
supported  by  an  aluminum-sealant-aluminum  sandwich  structure. 

Since  the  analysis  of  the  structure  of  Fig.  A-l(b)  Is  quite 
straightforward  for  SH  waves,  It  will  be  presented  here.  For  no  y-varlatlon 
and  x-propagatlon,  wc  require  the  stress  component  vanish  at  z  ■  tt/2. 

Also,  T4  and  the  particle  velocity  Vy  are  continuous  at  z  *  ± s/2 .  The 
boundary  conditions  are  separable,  hence  each  mode  satisfies  the  above 
boundary  conditions  Individually: 

Region  1 

Vy  *  An  cosk2n  (z  -  t/2) 
pV$2 

T4  '  '  TT  Ankzn  s1n  kzn  U  “  t/2) 


(1) 

(2) 
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Region  2 


vy  ■  Bn  “s  k'zn  z 


l  Symmetric 


p'Vj ? 

T4 - JT-k'zn  s1n  k'zn  2y 


(3) 

(4) 
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y 


T 
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B„  sin  k'  z 
n  zn 


cos 


■\ 

„  antisymmetric 


(5) 

(6) 


Region  3 

Vy  -  Cn  cos  kzn  (z  +  t/2)  (7) 

pvs2 

T4  "  ‘  TT  Cn  kzn  s1n  kzn  (z  '  t/2)  (8) 

and 

kx  "  '  kz2  f---  (9) 

Because  the  structure  Is  symmetrical ,  the  analysis  Is  straightforward.  By 
taking  the  ratio  T4/Vy  z  *  +  s/2  one  obtains  the  dispersion  relation  relating 
kx  to  u.  Cnee  this  relation  Is  known, then  Bn  and  Cn  are  given  simply  In  terms 
of  An  from  continuity  of  Vy  at  z  »  is/2  and  the  source  free  problem  Is 
solved. 
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Anti  symmetric 
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tt  -y 


(1  -  S/t) 


p'V, 


cot  IT 


f  S' 
an 

TT 


»  +  .0585  cot  In 


f  S' 
an 


The  limiting 


sn 


an 

n 


case  of  S  «  0  yields 

.  nV5 

~r 

(2n  +  1)  Vs 
- ^ - 

*  0,  1,  2,  3,  etc. 


(14) 

(15) 


(16) 

(17) 


which  corresponds  to  the  cutoff  frequencies  of  a  solid  aluminum  plate  of 
thickness  t.  As  S  Is  Increased  from  0,  but  remains  small  compared  to  t,  the 
cutoff  frequencies  are  shifted.  Analysis  shows  that  the  phase  velocity  of  the 

lowest,  order  symmetric  SH  mode,  which  has  no  cutoff  frequency,  Is  perturbed 

less  than  0.65%  for  t.  «■  250  and  S  <  10  mils.  This  mode  Is  the  analog  of  the 
n  *  0  mode  for  a  single  plate.  If  we  restrict  our  frequency  of  operation 
below  the  cutoff  frequency  for  the  next  mode,  fao,  corresponding  to  the  n  ■  1 
mode  of  a  simple  plate,  the  effects  of  sealant  variations  will  be 
negllgable.  Unfortunately,  from  (15),  one  finds  the  presence  of  a  soft 
material  such  as  polyethelene  or  sealant  drops  the  cutoff  frequency  to  the 
100-150  KHz  regime  for  S  on  the  order  of  2-5  mil  and  down  to  the  50-100  KHz 

regime  for  S  ~  10  mil.  We  know  EMATS  can  be  designed  In  the  range  but  we  must 

determine  If  resolution  becomes  a  limiting  factor  If  in  fact  one  needs  to 
operate  at  such  low  frequencies.  A  more  complete  characterization  of  the 
sealant  coupled  with  the  above  analysis  is  required. 
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